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摘要 : 天 体 视 位 置 是 进行 天 文 定 姿 的 必 备 数据 ， 目 前 视 位 置 计算 的 研究 多 基于 地 面 和 海 
上 ， 天 基 平 台 的 发 展 需 进一步 研究 。 介绍 了 星 上 视 位 置 转换 原理 ， 推 导 各 种 位 置 改正 公式 ， 基 
于 新 天 文 参 考 框架 建立 一 种 太阳 系 内 天 体 视 位 置 计 算 模型 ， 可 计算 任意 时 刻 天 体 的 地 心 视 位 置 
和 站 心 视 位 置 。 与 天 文 年 历 地 心 视 位 置 相 比 ， 类 木 行星 的 赤 经 差 值 在 10 ms 以 内 ， 杰 纬 差 值 在 
220 mas 以 内 ; 与 STK 仿真 的 站 心 视 位 置 相 比 ， 杰 经 差 值 在 0.15s 以 内 ， 杰 纬 差 值 在 2 以 内 。 
实验 表明 ， 该 模型 可 以 正确 计算 测 肯 太阳 系 内 天 体 的 视 位 置 ， 其 精度 满足 天 文 定 姿 的 需求 。 
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星 敏 感 器 利用 天 体 视 位 置 和 星空 图 像 天 量 信息 进行 天 文 定 姿 。 天 文 定 姿 利 用 恒星 、 日 月 以 及 行星 
= 等 自然 天 体 作为 导航 信 标 ， 确 定 航 天 咒 在 轨道 坐标 系 中 的 三 轴 姿 态 ， 天 文 定 姿 包 括 初始 模式 和 跟踪 模 
e 式 两 种 方式 '" 。 首 先 ， 星 敏感 器 对 视 场 内 天 体 成 像 ， 经 过 星 图 处 理 和 星 点 提取 等 技术 得 到 天 体 的 网 
- 像 矢量 信息 ， 对 星 图 上 的 星 点 数据 进行 星 图 匹配 ， 识 别 天 球 坐 标 系 中 对 应 的 天 体 ， 将 这 些 天 体 的 位 置 
r HR RIMI PAA D ROLE EE, KEEMIK MS BS ERR el fet I n] LU E RULES SB E, 
i 解 算 航 天 器 三 轴 姿 态 信息 。 因 此 ， 天 体 视 位 置 是 天 文 定 姿 的 必 备 数据 ， 为 提高 定 姿 精 度 ， 确 保 天 文 定 
姿 的 实时 性 要 求 ， 需 对 恒星 以 及 太阳 系 内 天 体 的 实时 视 位 置 计算 进行 研究 。 

目前 为 止 ， 国内 外 关于 天 体 视 位 置 计算 的 研究 主要 应 用 于 地 基 和 海 基 观测 。 文 [2] 基 于 DE405 历 
表 数据 ， 采 用 天 文 矢 量 处 理 方法 ， 提 出 了 一 种 高 精度 行星 视 位 置 计 算 的 数学 模型 ， 主 要 运用 于 舰 船 平 
台 。 文 [3] 归 纳 并 建立 了 常数 变易 法 和 坐标 摄 动 法 两 种 计算 行星 视 位 置 的 方法 ， 在 船舶 天 文 导航 系统 
中 使 用 。 文 [4-5] 介 绍 了 基于 地 面 观 测 的 新 型 天 文 测量 系统 中 行星 视 位 置 的 计算 方法 。 文 [6] 通 过 理 
论 分 析 得 出 一 种 行星 视 位 置 长 时 段 更 新 算法 ， 但 未 考虑 载体 运动 产生 的 影响 。 文 [7-8] 对 恒星 视 位 置 
的 计算 做 了 相关 研究 ， 对 于 太阳 系 内 天 体 视 位 置 的 计算 有 一 定 的 参考 作用 。 

天 文 定 姿 时 星 敏 感 顺 位 于 高 速 运动 的 航天 顺 上 ， 与 地 基 海 基 的 运用 情况 最 大 的 差别 就 是 测 站 本 刁 
产生 的 光 行 差 量 级 较 大 ， 不 可 忽略 (实验 结果 显示 影响 达 角 秒 量 级 )。 本 文 研究 了 新 天 文 参 考 框架 下 
天 文 定 姿 时 太阳 系 内 天 体 视 位 置 计算 ， 并 建立 了 一 种 实用 解 算 模型 。 与 2011 年 、2013 年 、2014 年 和 
2015 年 的 天 文 年 历 上 天 体 视 位 置 进行 比 对 ， 验 证 地 心 视 位 置 计算 的 精确 度 ， 利 用 STK 仿真 软件 模拟 
卫星 的 在 轨 运 行 ， 通 过 得 到 的 位 置 和 速度 信息 计算 天 体 的 序列 站 心 视 位 置 ， 与 STK 中 的 太阳 系 内 天 
体 视 位 置 数据 比 对 分 析 站 心 视 位 置 计算 的 准确 度 。 


1 太阳 系 内 天 体 视 位 置 转换 


1.1 新 天 文 参考 框架 
新 天 文 参考 框架 规定 了 新 的 坐标 系 和 岁差 章 动 模型 ， 天 体 视 位 置 计算 需 做 相应 的 改变 。 国 际 天 文 
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协会 (International Astronomical Union, IAU) 1991 年 引入 以 广义 相对 论 为 基础 的 国际 天 球 参 考 系 
( International Celestial Reference System, ICRS) ， 定 义 了 两 套 运 动 学 无 旋转 时 空 坐 标 系 太阳 系 质 心 天 球 
参考 系 ( Barycentric Celestial Reference System, BCRS) 和 地 心 天 球 参 考 系 (Geocentric Celestial Reference 
System, GCRS) , 2003 年 开始 采用 新 的 运动 瞬时 参考 框架 ,该 框架 同 动力 学 参考 架 没有 关系 ， 由 新 的 
岁差 章 动 模型 确定 。IAU2000 岁差 章 动 模型 描述 了 天 体 中 间 极 ( Celestial Intermediate Pole, CIP ) 在 地 心 
天 球 参考 系 中 的 运动 。 

(1)P03 岁差 模型 

天 文学 中 采用 过 几 种 岁差 模型 ，1986 年 以 前 一 直 使 用 纽 康 岁 差 系 统 ，1986 年 后 开始 采用 IAU1976 
岁差 模型 ，2000 年 后 采用 IAU2000 岁差 模型 ，2006 年 的 国际 天 文 协会 大 会 决定 2009 年 1 月 1 日 之 后 
采用 最 新 的 P03 模型 。P03 模型 的 岁差 量 如 (1) 式 
V, = 5038. 481507"T — 1. 079006"T^ — 0. 00114645"T^ + 0.000132851"T" - 9. 51" x 10^? 
w, = 84381. 406000" — 0. 025754"T + 0. 0512663"T? — 0. 00772503"T* 

-4.67 x 107"T' - 3.337" x 1077 
£, = 84381. 406000" — 46. 836769"T — 0. 0001831"T? + 0. 00200340"T* 

- 5.76" x 107 T* - 4.34" x 10* T? 
X, = 10. 556403"T - 2.381429"T? - 0. 00121197"T? + 0. 000170663"T" — 5. 60" x 10^? 
式 中 ,7 为 距离 历 元 J2000. 0 的 儒 略 世纪 数 。 

(2) IAU2000 章 动 模型 

1984 年 章 动 的 计算 采用 IAU1980 章 动 理论 ，2005 年 正式 采用 IAU2000 章 动 模型 。 其 中 ，IAU2000 章 
动 模型 分 为 2000A 和 2000B 模型 ，A 模型 精度 最 高 ， 优 于 0.2 mas， 它 包含 了 678 个 日 月 项 和 687 个 行星 
项 ， 总 共 1 365 项 ， 数 量 上 是 过 去 模型 的 十 多 倍 。 为 了 简化 计算 ,提高 解 算 效 率 ， 减少 计算 时 间 ， 对 A 
模型 进行 了 简化 ， 取 其 中 的 78 项 组 成 模型 B， 它 的 精度 在 1995-2050 年 之 间 可 以 保持 在 1 mas。 


678 


Ay -— £ (A; 十 4; 了) sinf, + (A;" T A:T) cosf. 
i=l 


(1) 


678 A 


Ae = 3, (B, + B/T)cosf, + (B," + B/"T)sinf, 
式 中 各 参数 详细 见 文 [11] 。 


i IA 岁 首 到 测 瞬 的 岁差 
1.2 视 位 置 转换 EPEE | EES 测 瞬 平 位 置 
DE/LE 历 表 里 包 含 了 太阳 系 内 所 有 天 体 的 ( 历 表 历 元 时 的 HUM ( 平 赤道 和 平 春分 点 
全 MR M i M ` 平 赤 道 坐 标 十 应 好 oy ZA 
基本 信息 ， 天 文 定 姿 时 太阳 系 内 天 体 的 测 瞬 站 心 a AHA E 


视 位 置 计算 需要 充分 考虑 各 种 因素 的 影响 。(1) 
天 体 的 自行 (proper motion): 从 岁 首 历 元 时 刻 到 
测 瞬时 刻 时 间 长 度 的 天 体 自行 ， 一 般 投影 到 星空 UE ELT 
中 自行 的 影响 小 于 1(/ 年 ); (2) EXE (parallax) : Sm 
两 个 不 同 地 点 看 到 的 物体 的 视 方 向 之 差 .， 视差 通 

常用 来 换算 天 体 之 间 的 距离 ， 它 的 量 级 小 于 1"; 
(3) 光 行 差 (aberration) : 由 观测 者 与 观测 对 象 之 


章 动 改正 
(EOP 参 数 ) 


视差 改正 
光 引 力 偏 折 改 正 
光 行 差 改正 


间 的 相对 运动 和 光速 有 限 引 起 的 位 置 视觉 差异 ，。 | “站 心 视 位 置 “| exert | 。 地 心 视 位 置 
MBA ` NP 全 jé TESA 真 赤 道 和 真 春 分 点 真 赤道 和 真 春 分 点 

在 地 面 进行 天 体 的 观测 时 光 行 差 的 影响 可 以 达到 | Gassen S T 

21"， 而 航天 器 在 轨道 上 的 运行 速度 对 视 位 置 的 u 

计算 影响 更 大 ， 是 研究 航天 器 上 恒星 视 位 置 计算 5 NM 

与 地 面 观测 之 间 的 主要 差别 (4) 还 有 光 引 力 偏 Fig. 1 Relationship between various apparent positions 

puoi A y ; ki J 


折 、 岁 差 和 章 动 等 。 图 1 给 出 了 天 文 定 姿 中 天 体 各 种 位 置 之 间 的 转换 关系 。 
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2 视 位 置 计算 模型 


同 恒星 的 视 位 置 计算 相 比 ， 太 阳 系 内 天 体 的 视 位 置 计算 最 大 的 区 别 在 于 需要 考虑 由 光 行 差 引 起 的 
位 置 影响 。 因 为 测 站 同 被 观测 天 体 距离 近 ， 相 对 运动 速度 大 ， 故 需要 考虑 周 日 视差 和 周 日 光 行 差 改 
正 。 具 体 计算 步骤 如 下 : 

(1) 计 算 观 测 瞬 间 的 地 球 时 TT 对 应 的 侨 略 日 J 和 TDB 时 间 

鉴于 早期 使 用 的 历 书 时 存在 不 连续 和 不 一 致 等 问题 以 及 原子 频 标 技术 的 高 速 发 展 ，1976 年 起 采用 力 
学 时 以 广义 相对 论 作为 时 间 工 作 的 理论 基础 。 目 前 采用 两 种 力学 时 ， 太 阳 系 质心 参考 系 用 质心 力学 时 
7DB， 地 心 视 位 置 的 解 算 用 地 球 时 TT, TT 是 建立 在 国际 原子 时 TAI 的 基础 上 ， 两 者 关系 如 式 : 

TT = TAI + 32.184 = UTC + LS + 32.184, (3) 

LS 为 跳 秒 ， 可 根据 IERS 上 面 发 布 的 公报 B 参数 数据 确定 跳 秒 的 值 。 

(2) 通 过 DE/LE 历 表 内 捕 得 到 J2000 历 元 下 太阳 系 内 各 天 体质 心 的 位 置 拓 量 rT) 和 地 球 在 TDB 
时 刻 对 应 的 质心 天 球 参考 系 中 的 质心 位 置 r, MEE v., v, 是 地 球 绕 太 阳 的 公转 速度 ， 单 位 为 AU/ 日 。 

(3) 计算 航天 器 TT 时刻 在 真 赤道 真 春分 点 坐标 体系 下 相对 于 地 心 的 位 置 r, RUE BE v., r, 以 km 
ARE, v, 以 km/s 为 单位 ， 从 而 得 到 航天 器 在 质心 天 球 参考 系 中 的 位 置 > 和 速度 v, 


r,-r,tr,., 


v, =v, ty. (4) 
> (4) 计 算 太 阳 系 天 体 的 光 行 时 。 光 行 时 的 近似 值 > 为 
o r= |r(2 -r, le. (5) 
- 根据 DE405 历 表 获 取 7-7 时 太阳 系 内 天 体 的 质心 位 置 矢量 r(T-7) ， 则 光 行 时 改进 值 为 
ham 7’=|r(T-7) espe (6) 
Jin 依 此 循环 交代 ， 直 到 光 行 时 满足 精度 要 求 。 
€ (5) YE T-r 时 太阳 系 内 天 体 相 对 于 了 时 航天 器 的 坐标 方向 的 单位 矢量 己 ; 
r,( T-T7)-r, 
P nim 
p (7) 
p=|r,(T-7) -r,l 
(6) 进 行 光线 引力 弯曲 改正 、 周 年 光 行 差 改 正 、 岁 差 和 章 动 改正 : 
jg c e Pe. 
s clr 9d X r, p 
1 P-V, 1 (8) 
P, =P, 和 er" ):Px(PxV, so ates (PxV), 
P,=N-P-P,. 
(7) 将 位 置 矢量 转换 为 球面 坐标 : 
a, = tan (y/x,) ， 
(9) 
8, 2 tan  (z,/ Mx3 + yi). 
(8) 加 入 周 日 视差 改正 : 
a, = a, — Pucospsec6sint， 
1 "T" (10) 
ô, =ô, — singcosó, + cosqQsinÓ,cost + ru — a ) sinó,sint . 
P, 为 行星 的 周 日 地 平视 差 ， 可 根据 历 表 获 取 ; e 为 测 站 天 文 纬度 ; t 为 天 体 的 时 角 。 
(9) 加 入 周 日 光 行 差 改 正 : 
a = a, — (^. 021coseoostsecó, , 
(11) 


8 =ô, — 0. 32cosgsinó,sint . 
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利用 VS2010 平台 对 太阳 系 内 天 体 视 位 置 计算 进行 了 编译 ， 历 表 采 用 美国 喷气 推进 实验 室 提供 的 
DE405 历 表 数 据 ， 岁 差 章 动 模 型 采用 最 新 的 IAU2000 岁差 章 动 模型 。 天 体 视 位 置 解 算 程 序 界面 如 图 2。 
在 参数 输入 处 输入 测 瞬 历 元 信息 、 测 站 信息 和 对 应 的 地 球 定向 参数 (Earth Orientation Parameter, 


EOP), ， 选 择 对 应 的 解 算 类 型 。 表 1 地 心 视 位 置 计算 差 值 
Table 1 The difference value of geocentric 
3 Ki ij 计算 和 精度 分 析 apparent position 
天 体 名 称 杰 经 差 值 /s 赤 纬 差 值 /” 

3.1 地 心 视 位 置 计算 Mercury -0. 1004 -1.3312 

分 别 采 用 2011、2013、2014 和 2015 年 4 个 Venus -0.0217 0 099 0 
年 份 的 天 文 年 历 验证 程序 的 准确 性 ， 观 测 时 间 取 VR ovis Nt 
每 年 的 2 月 3 日 力学 时 0 时 ， 测 站 统一 设置 在 地 心 rel NOM i ge. 
处 。 限 于 篇 幅 ， 表 1 仅 列 出 2014 年 的 比 对 数据 。 S M d 

从 表 1 可 以 看 出 ， 离 航天 器 越 远 的 行星 ， 计 MAR boss Pos 
算 的 精度 越 高 而 太阳 、 水 星 、 金 星 和 火星 等 这 Neptune —0. 001 7 -0. 0563 
些 在 内 行星 范围 内 的 类 地 行星 其 计算 精度 比较 an epe eim 


低 ， 且 离 航天 器 位 置 越 近 精 度 越 低 ， 类 地 行星 的 
赤 经 误差 在 0.35 s UN, WRZE 2 以内。 类 木 行星 的 赤 经 误差 在 10 ms 以 内 ， 赤 纬 误 差 在 220 mas 
以 内 。 图 3 给 出 了 4 年 计算 结果 同 天 文 年 历 比 对 的 差 值 ， 天 体 序 列 顺序 同 表 1。 


gxsamem A E """Ww* i8 x 
参数 办 
GEARS 历 表 信 息 
日 其 
2015 年 02 月 09 日 DE405 
3 时 12 分 4534 W 
测 站 信息 选择 解 算 类 型 
x 6321 y -2100 I 231 
地 心 视 位 置 
wal vw X3  v 2 
大 气压 : 1001 温度 14 
e 测 站 视 位 置 
EOS 
Wb: 35 VT1-UIC 0.4412 
开始 解 算 
= 行星 
图 2 视 位 置 解 算 界 面 图 3 地 心 视 位 置 解 算 差 值 
Fig.2 The interface of apparent position calculation Fig.3 The difference value of geocentric apparent position 


3.2 ”站 心 视 位 置 计算 

利用 STK 卫星 轨道 仿真 软件 ， 创 建 一 颗 600 km 高 度 的 人 造 卫 星 ， 考 虑 J2 摄 动 确定 定 轨 数 据 ， 轨 
道 倾角 i= 60.074 deg， 升 交点 赤 经 2=298.074 deg， 偏 心率 设 为 0， 了 卫星 的 姿态 信息 采用 经 典 模式 。 
模拟 2015 年 2 月 7 日 到 8 日 两 天 的 轨道 数据 ， 时 间 间 隔 取 5 min。 通 过 STK 仿真 得 到 577 组 卫星 的 位 


置 和 速度 信息 ， 代 和 人 天 体 视 位 置 程序 计算 站 心 视 X2 站 心 视 位 置 计 算 指标 
位 置 。 计 算 J2000 坐标 系 下 地 球 、 月 球 和 太阳 的 Table 2 The counting indicator of observer-fixed 
视 位 置 ， 得 到 赤 经 赤 纬 差 值 如 图 4。 apparent position 

由 图 4 可 看 出 赤 经 赤 纬 差 值 较 小 ， 赤 经 差 值 Aa Aa A8 A8 
在 0.15s 以内， 赤 纬 差 值 在 2 以 内 ， 从 侧面 验 E ast. 均 方 差 /。 。 均值” 均 方差 / 
证 了 站 心 视 位 置 计 算 的 正确 性 ，STK 所 提供 的 视 Earth -0.0012 0.0709 -0.0133 1.0134 
位 置 也 不 是 最 精确 模型 解 算 所 得 ， 因 此 该 结果 不 Moon 0.0325 0.0687 -0.1926 1.0234 


能 作为 天 体 视 位 置 计算 程序 的 精度 。 表 2 给 出 了 Sw 0.0019 0.0693 0.0199 1.0577 
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地 球 、 月 球 和 太阳 站 心 视 位 置 解 算 的 赤 经 赤 纬 均值 和 均 方差 .可 以 看 出 ， 地 球 的 视 位 置 计算 精度 最 
低 ， 因 为 它 距 离 航天 器 最 近 。 太 阳 的 赤 经 差 值 均值 1. 9 ms ， 赤 纬 均值 19. 9 mas， 在 白天 观 星 困 难 的 情 
况 下 ， 太 阳 是 最 为 稳定 的 信息 源 ， 高 精度 太阳 视 位 置 计算 对 天 文 定 姿 具 有 很 好 的 帮助 作用 。 
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图 4 ”站 心 视 位 置 解 算 差 值 
Fig.4 The difference value of observer-fixed apparent position 
4 结 论 


本 文 研究 了 天 文 定 姿 时 太阳 系 内 天 体 视 位 置 的 计算 方法 ， 建 立 一 种 实用 的 视 位 置 算 法 模型 ， 该 模 
型 可 以 计算 天 体 的 地 心 视 位 置 和 站 心 视 位 置 。 与 2011、2013 2014 和 2015 年 的 天 文 年 历 比较 发 现 ， 
类 地 行星 的 赤 经 误差 在 0.35 s 以 内 ， 赤 纬 误差 在 2 以 内 ， 类 木 行星 的 赤 经 误差 在 10 ms 以 内 ， 赤 纬 误 
差 在 220 mas 以 内 ， 验 证 了 地 心 视 位 置 计算 的 准确 性 。 利 用 STK 仿真 软件 模拟 在 轨 卫 星 两 天 的 轨道 数 
据 ， 计 算 地 球 、 月 球 和 太阳 的 站 心 视 位 置 ， 与 STK 得 到 的 J2000 坐标 系 下 的 视 位 置 对 比 ， 赤 经 差 值 在 
0.15 s 以 内 ， 赤 纬 差 值 在 2” 以 内 ， 太 阳 的 赤 经 差 值 均值 1. 9 ms， 赤 纬 均值 19. 9 mas ， 验 证 了 站 心 视 位 
置 计算 的 正确 性 。 由 于 在 天 文 单位 中 ， 赤 经 的 使 用 单位 是 时 秒 ， 而 赤 纬 单位 是 角 秒 ， 时 秒 与 角 秒 之 间 
相差 15 倍 ， 若 进行 单位 转换 之 后 ， 可 以 看 出 视 位 置 计算 的 赤 经 、 赤 纬 误差 大 小 的 量 级 差不多 。 实 验 
表明 ， 该 计算 模型 可 以 正确 解 算 测 瞬 太 阳 系 内 天 体 的 视 位 置 ， 精 度 满足 天 文 定 姿 的 需求 。 
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Apparent Position Calculation of Celestial Body in the Solar 
System Used for the Celestial Attitude Determination 
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Abstract: Celestial apparent position is the key data for celestial attitude determination. At present, the 
study of apparent position calculation mainly focuses on ground and ocean; development in space-based platform 
should take a further step. A conversion principle of apparent position in the spacecraft and various position 
correction formula is introduced in this thesis; based on the new astronomical reference frame, a model is 
established to calculate the apparent position of celestial bodies in the solar system, through which instantaneous 
geocentric virtual position and observer-fixed apparent position can be calculated. Compared with geocentric 
virtual position of astronomical almanac, the difference value of Jovian planet’ right ascension can be controlled 
within 10ms and declination difference under 220mas; Compared with observer-fixed apparent position of 
satellite tool kit, the difference value of right ascension can be kept under 0. 15s and the declination difference 
under 2". Experiments show that the model can correctly compute instantaneous apparent position of celestial 
body in the solar system and such precision meets the requirements of celestial attitude determination. 


Key words: Apparent position calculation; Celestial attitude determination; Astronomical almanac ; System Tool Kit 


